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In this paper, we continue our development of algorithms used for topological network discovery. We 
present native P system versions of two fundamental problems in graph theory: finding the maximum 
number of edge- and node-disjoint paths between a source node and target node. We start from the 
standard depth-first-search maximum flow algorithms, but our approach is totally distributed, when 
initially no structural information is available and each P system cell has to even learn its immediate 
neighbors. For the node-disjoint version, our P system rules are designed to enforce node weight 
capacities (of one), in addition to edge capacities (of one), which are not readily available in the 
standard network flow algorithms. 
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1 Introduction 

Inspired by the structure and interaction of living cells, P systems provides a distributed computational 
model, as introduced by G. Paun in 1998 [12]. The model was initially based on transition rules, but was 
later expanded into a large family of related models, such as tissue and neural P systems (nP systems) 
t71 [T3l and hyperdag P systems (hP systems) [ 10]. Essentially, all versions of P systems have a structure 
consisting of cell-like membranes and a set of rules that govern their evolution over time. A large 
variety of rules have been used to describe the operational behavior of P systems, the main ones being: 
multiset rewriting rules, communication rules and membrane handling rules. Transition P systems and 
nP systems use multiset rewriting rules, P systems with symport/antiport operate by communicating 
immutable objects, P systems with active membranes combine all three type rules. For a comprehensive 
overview and more details, we refer the reader to |[T3Tl . 

Earlier in [3], we have proposed an extensible framework called P modules, to assist the programma- 
bility of P systems. P modules enable the modular composition of complex P systems and also embrace 
the essential features of a variety of P systems. In this paper, we will use a restricted subset of this uni- 
fying model, called simple P modules, (subset equivalent to neural P systems [7 ]), to develop algorithms 
for finding the maximum number of edge- and node-disjoint paths between two cells in a fairly large 
class of P systems, where duplex communication channels exist between neighboring cells. We assume 
that the digraph structure of the simple P module is completely unknown (even the local neighboring 
cells must be discovered O) and we need to, via a distributed process, optimally create local routing 
tables between a given source and target cell. 

There are endless natural applications that need to find alternative routes between two points, from 
learning strategies to neural or vascular remodeling after a stroke. In this paper, we focus on a related 
but highly idealized goal, how to compute a maximum cardinality set of edge- and node-disjoint paths 
between two arbitrary nodes in a given digraph. 
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One obvious application related to networks is to find the best bandwidth utilization for routing of 
information between a source and target lfl~5ll . For instance, streaming of applications for multi-core 
computations uses edge-disjoint paths routing for task decomposition and inter-task communications 
ifTTl . In fact, classical solutions are based on a network flow approach such as given in (21 |H, or on 
Menger's Theorem, an old, but very useful, result, cited below. 

Theorem 1 (Menger (8l). Let D = (V,A) be a digraph and let s,t EV. Then the maximum number of 
node-disjoint s—t paths is equal to the minimum size of an s—t disconnecting node set. 

Another application is to find a maximum matching (or pairing) between two compatible sets such 
as the marriage arrangement problem or assigning workers to jobs. 

Our third application (and a motivating problem for the authors) is the Byzantine Agreement problem 
EllH, in the case of non-complete graphs. The standard solution (also based on Menger's Theorem) 
allows for k faulty nodes (within a set of nodes of order at least 3k + 1) if and only if there are at 
least 2k + 1 node-disjoint paths between each pair of nodes, to ensure that a distributed consensus can 
occur |6]. 

Briefly, the paper is organized as follows. In the next section, we give a formal definition of sim- 
ple P modules, to give a unified platform for developing our P systems algorithms. Next, in Section|3]we 
summarize the standard network flow approaches for finding edge- and node-disjoint paths in digraphs 
and we discuss optimizations and alternative strategies which are more appropriate for P systems. In 
Section |4| we discuss three possible relations between the structural digraph underlying a simple P mod- 
ule and the search digraph used for determining paths. Section [5] details breadth-first-search rules used 
to determine the local cell topologies, i.e. all cell neighborhoods; this is a common preliminary phase for 
both the edge- and node-disjoint path implementations. The next two sections detail depth-first-search 
rules for the edge-disjoint case (in Section [6]) and for the node-disjoint case (in Section [7]). Finally, in 
Section [8] we end with conclusions and some open problems. 

2 Preliminary 

We assume that the reader is familiar with the basic terminology and notations: functions, relations, 
graphs, edges, nodes (vertices), directed graphs, arcs, paths, directed acyclic graphs (dags), trees, alpha- 
bets, strings and multisets ifTTTl . We now introduce simple P modules, as a unified model for representing 
several types of P systems. Simple P modules are a simplified variety of the full P modules, which omit 
the extensibility features and use duplex communication channels only. With these restrictions, although 
their formal definitions are different, simple P modules are essentially equivalent to neural P systems 0. 
For the full definition of P modules and further details on recursive modular compositions, the reader is 
referred to . 

Definition 2 (simple P module). A simple P module is a system IT = (0,K, 8), where: 

1. O is a finite non-empty alphabet of objects; 

2. K = {oi , 02, . . . , <7„} is a finite set of (internal) cells; 

3. 8 is a binary relation on K, without reflexive or symmetric arcs, which represents a set of parent- 
child structural arcs between existing cells, with duplex communication capability. 

Each cell, a, E K, has the initial form a, = (Qi,so,wo,Ri) and general form a, = (Qi,s,w,Ri), where: 

• Qj is a finite set of states; 



M.J. Dinneen, Y.-B. Kim & R. Nicolescu 



123 



• s o £ Qi is the initial state; s G £?, is the current state; 

• wo € 0* is the initial multiset of objects; w G O* is the current multiset of objects; 

• is a finite ordered set of multiset rewriting rates of the general form: s x — s' x' (u)^ , where 
s,J G Q, x,x' eO*,ue O*, a G {min,max}, /3 G {t>4->$}. /G {one, spread, repl}U A". If u = X, 
denoting the empty string of objects, this rule can be abbreviated as s x — > a s' x'. The application 
of a rule takes two sub-steps, after which the cell's current state s and multi-set of objects x is 
replaced by s' and x', respectively, while u is a message which is sent as specified by the transfer 
operator j3 r . 

The rules given by the ordered set(s) Rj are applied in the weak priority order (I4ll . For a cell 
Ci = (Qht,w,Rj), a rule s x — > a s' x' (u)^ G Rt is applicable ift = s and x C w. Additionally, if s x — > a 
/ (m)^ is the first applicable rule, then each subsequent applicable rule's target state (i.e. state indicated 
in the right-hand side) must be s'. The semantics of the rules and the meaning of operators a, j8, 7 are 
now described. 

For convenience, we will often identify a cell a, with its index (or cell ID) i, when the context of 
the variable i is clear. We accept that cell IDs appear as objects or indices of complex objects. Also, we 
accept custom cell ID rules, which distinguish the cell ID of the current cell from other cell IDs. For 
example, the rule 0. 1 of Section[5] given as "s$ gi — s- min so" f° r ce U &u appears as "so g\ — s-min V in cell 
G{ and as "s gi ~W s " in cell a 2 . 

The rewriting operator a = max indicates that an applicable rewriting rule of /?, is applied as many 
times as possible, while the operator a = min requires a rule of /?, is applied only once. The communica- 
tion structure is based on the underlying digraph structure. In this paper, and we will only use the j8 = £ 
and 7 G {repl} UK transfer operators. With reference to cell a,-, a rewriting rule using («)^ repl indicates 
that the multiset u is replicated and sent to all neighboring cells (parents and children), i.e. to all cells in 
8(f) U 8 (i). Assuming that cell Oj is a parent or a child of cell a,-, i.e. j G 8(i) U a rewriting 

rule using (h)^ indicates that the multiset u is specifically sent cell Oj. Otherwise, if j ^ 8(i) U 8~ l (i), 
the rule is still applied, but the message u is silently discarded. The other non-deterministic transfer 
operators (e.g., one, spread, y, |) are just mentioned here for completeness, without details, and are not 
used in this paper. For details, the interested reader is referred to (31 • 

Remark 3. This definition of simple P module subsumes several earlier definitions of P systems, hP sys- 
tems and nP systems. If 8 is a tree, then II is essentially a tree-based P system (which can also be 
interpreted as a cell-like P system). If 8 is a dag, then II is essentially an hP system. If 8 is a digraph, 
then II is essentially an nP system. 

3 Disjoint paths in digraphs 

We now briefly describe the basic edge- and node-disjoint paths algorithms, based on network flow, 
particularized for unweighted edges (i.e. all edge capacities are one), see Ford and Fulkerson [5\. Our 
presentation will largely follow the standard approach, but also propose a couple of customizations and 
optimizations, specifically targeted for running on highly distributed and parallel computing models, 
such as P systems. 

We are given a digraph G = (V,E) and two nodes, a source node, s G V, and a target node, t G V. We 
consider the following two optimization problems: (1) find a maximum cardinality set of edge-disjoint 
paths from s to t; and (2) find a maximum cardinality set of node-disjoint paths from s to t. Obviously, 
any set of node-disjoint paths is also edge-disjoint, but the converse is not true. For example: 
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• Figure[T](a) shows a maximum cardinality set of edge-disjoint paths for a digraph G, which is also 
a maximum cardinality set of node-disjoint paths, 

• Figure [TJ (b) shows two maximum cardinality sets of edge-disjoint paths for the same digraph G, 
which are not node-disjoint. 

• Figure [2] shows a digraph where the maximum number of edge-disjoint paths is greater than the 
maximum number of node-disjoint paths. 




Figure 1: For this digraph, the maximum number of edge-disjoint paths from s to t, which is 2, can be 
achieved in three ways: (a) paths set {s.w.x.t, s.y.z-t}; (b) either of the following two paths sets: {s.w.x.t, 
s.y.x.z.t}, {s.w.x.z.t, s.y.x.t}. Paths shown in (a) are also node-disjoint, but paths shown in (b) are not. 




Figure 2: For this digraph, the maximum number of edge-disjoint paths from s to t (2) is greater than the 
maximum number of node-disjoint paths (1). 



3.1 Edge disjoint paths in digraphs 

In both edge- and node-disjoint cases, the basic algorithms work by repeatedly searching paths, called 
augmenting paths, in an auxiliary structure, called residual network or residual digraph. We will first 
focus more on edge-disjoint paths, because the node-disjoint paths can be considered as an edge-disjoint 
paths problem, with additional constraints. 

For the following "network flow" definition for digraphs with non-weighted arcs, we say that an arc 
(u,v) is in a set of paths P, denoted by the slightly abused notation (u,v) G P, if there exists a path % G P 
that uses arc (w, v). 

Definition 4. Consider a digraph G = (V,E), two nodes s and t, {s,t} C V, and a set P of edge-disjoint 
paths from s to t. Nodes in P are called flow-nodes and arcs in P are called flow-arcs. 

Given path % G P, each flow-arc (u, v) G % has a natural incoming and outgoing direction-the flow is 
from the source to the target; with respect to 7t, u is the flow-predecessor of v and v is the flow-successor 
of u. 
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The residual digraph is the digraph R = (V,E'), where the arcs in P are reversed, or, more formally, 
E' = (E \{(w,v) I (k,v) G P}) U {(v,u) I (k,v) G P}. Any path from 5 to t in /? is called an augmenting 
path. 

Given augmenting path a, each flow-arc (u, v) G a has also a natural incoming and outgoing direction- 
the flow is from the source to the target; with respect to a, u is the search-predecessor of v and v is the 
search-successor of m. 

Tfacf 5. Augmenting paths can be used to construct a larger set of edge-disjoint paths. More precisely, 
consider a digraph G and two nodes s and t. A set Pk of k edge-disjoint paths from s to t and an 
augmenting path a from s to ? can be used together to construct a set P^ +1 of k + 1 edge-disjoint paths. 
First, paths in {a} UP& are fragmented, by removing "conflicting" arcs, i.e. arcs that appear in QUQ, 
where Q = PHa (where ~ indicates arc reversal). Then, new paths are created by concatenating resulting 
fragments. For the formal definition of this construction, we refer the reader to Ford and Fulkerson Q. 
Note that including a reversed arc in an augmenting path is known as^ow pushback operation. 

This construction is illustrated in Figure [3] Figure [3] (a) illustrates a digraph G and a set Pi of 
edge-disjoint paths from s to t, currently the singleton {?To}, where 7To = s.y.x.t. Figure |3](b) shows its 
associated residual digraph R (note the arcs reversal). Figure [3] (c) shows an augmenting path a in R, 
a = s.w.x.y.z.t. Figure[3](d) shows the extended set P2 (after removing arcs (x,y) and (y,x)), consisting 
of two edge-disjoint paths from s to t, %\ = s.w.x.t, %2 = s.y.z-t. Figure [4] shows a similar scenario, 
where another augmenting path is found. Note that the two paths illustrated in Figure [3] (d) form both 
a maximum edge-disjoint set and a maximum node-disjoint set; however, the two paths sets shown in 
Figure [4] (d) form two other maximum edge-disjoint path sets, but none of them is node-disjoint. 




(a) (b) (c) (d) 

Figure 3: A residual digraph and an augmenting path: (a) a digraph G and one (edge-disjoint) path 7To 
from s to t (indicated by bold arrows), (b) the residual digraph Ro associated to digraph G and path 71q. 
(c) an augmenting path a in Rq (indicated by hollow arrows), (d) two new edge-disjoint paths 7t\ and %2, 
reconstructed from 71q and a (both indicated by bold arrows). 




Figure 4: The residual digraph of Figure [3] with another augmenting path and two new paths sets, 
{s.w.x.t,s.y.x.z.t}, {s.w.x.z.t, s.y.x.t}, which are edge-disjoint but not node-disjoint. 
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The pseudo-code of Algorifhm[6]effectively finds the maximum number (and a representative set) of 
edge-disjoint paths from s and t. 

Algorithm 6 (Basic edge-disjoint paths algorithm). 

1 Input : a digraph G = (V,E) and two nodes s G V, t G V 

2 k = (the stage counter) 

3 Po = Q (the current set of edge-disjoint paths) 

4 Rq = G (the current residual digraph) 

5 loop 

6 a = an augmenting path in R^, from s to t, if any (this is a search operation) 

7 if a = null then break 

8 k = k + 1 (next stage) 

9 Pjt = the larger paths set constructed using P,t-i and a (as indicated in Fact[5j) 

10 Rk = the residual digraph of G and P k 

11 end loop 

12 Output : & and P^, i.e. the maximum number of edge-disjoint paths and a representative set 

Typically, the internal implementation of search at step 6 alternates between a forward mode in 
the residual digraph, which tries to extend a partial augmenting path, and a backwards backtrack mode 
in the residual digraph, which retreats from an unsuccessful attempt, looking for other ways to move 
forward. The internal implementation of step 9 (i.e. Fact[5]> walks backwards in the residual digraph, as a 
consolidation phase, which recombines the newly found augmenting path with the existing edge-disjoint 
paths. 

This algorithm runs in k+ 1 stages, i.e. in up to outdegree(^) + 1 stages, if we count the number 
of times it looks for an augmenting paths, and terminates when a new augmenting path is not found. 
The actual procedure used (in step 6) to find the augmenting path separates two families of algorithms: 
(1) algorithms from the Ford-Fulkerson family use a depth-first-search (DFS); (2) algorithms from the 
Edmonds-Karp family use a breadth-first-search (BFS). As usual, both DFS and BFS use "bread crumb" 
objects, as markers, to avoid cycles; at the end of each stage, these markers are cleaned, to start again 
with a fresh context. In this paper, we develop P algorithms from the Ford-Fulkerson family, i.e. using 
DFS. 

3.2 Node disjoint paths in digraphs 

The edge-disjoint version can be also used to find node-disjoint paths. The textbook solution for the 
node-disjoint problem is usually achieved by a simple procedure which transforms the original digraph 
in such a way that, on the transformed digraph, the edge-disjoint problem is identical to the node-disjoint 
problem of the original digraph. Essentially, this procedure globally replaces every node v, other than 
s and t, with two nodes, an entry node vi and an exit node v 2 , connected by a single arc (vi,V2). More 
formally, the new digraph G' = (V',E') has V' = {s,t} U {vi,v 2 | v € V \ {s,t}, E' = {(vi,v 2 ) | v G 
V \ {s,t}} L){(u2, vi) I (u, v) G E}, where, for convenience, we assume that si =S2 = s and t\ = ? 2 = t are 
aliases. This standard node-splitting technique is illustrated in Figure [5] It is straightforward to see that 
the newly introduced arcs (wi,w 2 ), (xi,x 2 ), (^1,^2) and (zi,z 2 ), constrain any edge-disjoint solution to 
be also node-disjoint. 

However, in our case, since each node is identified with a P systems cell, we cannot solve the node- 
disjoint paths problem using the standard node-splitting technique. We propose two non-standard search 
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Figure 5: The node splitting technique. 



rules, which together limit the out-flow capacity of each v G V \ {s,t} to one, by simulating the node- 
splitting technique, without actually splitting the nodes. We believe that our rules can be used in other 
distributed network models where the standard node-splitting technique is not applicable. These rules 
are illustrated by the scenario presented in Figure [6j where we assume that we have already determined 
a first flow -path, s.x.y.z.t, and we are now trying to build a new augmenting path. 

1 . Consider the case when the augmenting path, consisting of s, tries flow-node y via the non-flow arc 
(s,y). We cannot continue with the existing non-flow arc (y,t) (as the edge-disjoint version would 
do), because this will exceed node y's capacity, which is one already. Therefore, we continue the 
search with just the reversed flow-arc (y,x). Note that, in the underlying node-splitting scenario, 
we are only visiting the entry node y\, but not its exit pair y2- 

2. Consider now the case when the augmenting path, extended now to s.y.x.z, tries again the flow- 
node y, via the reversed flow-arc (z,y). It may appear that we are breaking the traditional search 
rules, by re-visiting the already visited node y. However, there is no infringement in the underlying 
node-splitting scenario, where we are now trying the not-yet-visited exit node j2 (to extend the 
underlying augmenting path s.y\.X2.Zi). From y, we continue with any available non-flow arc, 
if any, otherwise, we backtrack. In our example, we continue with arc (y,t). We obtain a new 
augmenting "path", s.y.x.z,.y.t (corresponding to the underlying augmenting path s.y\.X2.z.\-y2-t)- 
We further recombine it with the already existing flow-path s.x .y.z.t, and we finally obtain the two 
possible node-disjoint paths, s.x.z.t and s.y.t. 




Figure 6: Node-disjoint paths, (a) non-standard search: flow path s.x.y.z.t and augmenting "path" 
s.y.x.z.y.t. (b) node-splitting: flow path s.x\.X2.y\.y2-Z\-Z2-t and augmenting path s.yi.X2-Z\.y2-t- 

The following theorem is now straightforward: 

Theorem 7. If the augmented path search in step 6 of Algorithm^is modified as indicated above, the 
algorithm will terminate with a restricted subset of edge-disjoint paths, forming a maximum cardinal 
subset of node-disjoint paths. 
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3.3 Pointer management 

With respect to the implementation, the edge-disjoint version provides its own additional challenge, 
not present in the node-disjoint version. In the node-disjoint version, each flow-node needs only one 
pointer to its flow-predecessor and another to its flow-successor. However, in the edge-disjoint version, 
a flow-node can have k flow -predecessors and k flow-successors, with k > 1, where each combination is 
possible, giving rise to k\ different edge-disjoint paths sets, each of size k, passing through this node. A 
naive approach would require recording full details of all k\ possible size-& paths sets, or, at least, full 
details for one of them. 

In our simplified approach, we do not keep full path details; instead, a node needs only two size 
k lists: its flow-predecessors list and its flow-successors list. Using this information, any of the actual 
k\ paths sets can be recreated on the fly, by properly matching flow-predecessors with flow-successors. 
As an example, consider node x of Figure |4] (d), which has two flow-predecessors, w and y, and two 
flow-successors, t and z; thus w is part of four distinct paths. Node w needs only two size-fc lists: its 
flow -predecessors list, and its flow-successors list {z,t}- 

3.4 Possible optimization 

We propose a potential speed-up for Algorithm [6j We restrict this discussion to edge-disjoint paths 
and standard DFS; however, the discussion can be generalized to more general flows and other search 
patterns. Consider V s = En ({s} x V) = {(s, vi), (s, V2), . ■ ■ (s, VrfJ}, where d s = outdegree(s). Step 6 
systematically tries all arcs in V s . Without loss of generality, we assume that step 6 always tries arcs in the 
order (j,vi), (s,V2), ■ ■ ■ (s,v<i s ), stopping upon the first arc which is identified as starting an augmenting 
path a, say (s,v r ), where r E For brevity, we will indicate this by saying that arc (s,v r ) is the first 

one that succeeds (and arcs (s,vi), (s,V2), ■ ■ ■ (s, v r -\)fail). 

Consider a complete run of Algorithm [6j Assume that this algorithm finds k augmenting paths and 
then stops. Assume that stage j G [l,k], finds a new augmenting path ctj, which starts with arc (s,v, ; ), 
i.e. arc (s,v,v) is the first one that succeeds at stage j. A direct implementation of Algorithm |6] seems to 
require that step 6 starts a completely new search for each stage, restarting from (s, vi) and retrying arcs 
that have been previously tried (whether they failed or succeeded). 

However, this is not necessary. Theorem[8]indicates that stage j + 1 does not need to retry the nodes 
that have already been considered (whether they failed or succeeded). Specifically, the indices indicating 
the successful arcs in V s are ordered by stage number, z'i < ii < ••• < 4, and, at stage j € [l,k], we can 
start the new search directly from arc (s,v,v (where z'q = 0). 

Theorem 8. Step 6 of Algorithm^can start directly with the arc in V s which follows the previous stage's 
succeeding arc in V s . 

Proof. At stage j G [l,k], after finding augmenting path ctj, the algorithm fragments a, together with 
the previous set of edge-disjoint paths, Pj, deletes some arcs and reassembles a new and larger set of 
edge-disjoint-paths, Pj + i, such that |^+i| = \Pi\ + 1. 

(1) Consider arc (s, v;.), the starting arc of path ctj. We first show that a successful arc, such as (j,v;.), 
need not be tried again by step 6, for two arguments: 

(1.1) Arc (s,Vi ) cannot start another augmenting path, because all following residual digraphs will 
only contain its reversal (v,v,s), never its direct form (s, v,- ; ). 

(1.2) Arc (s,Vj ) cannot be revisited as part of another augmenting path. No augmenting path con- 
tains arcs from V x {s}, because the search (DFS or BFS) avoids already visited nodes (marked with 
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"pebbles"), and s is always the starting point and thus the first node marked. Therefore, once successful 
arc (s,v ( - ) is never deleted by "flow pushback operations". 

From (1.1) and (1.2), we conclude that, once successful, an arc in V s will always be the starting arc 
of an edge-disjoint path and need not be tried again by step 6. 

(2) We next show that, once failed, an arc in V s will always fail. This part of the proof is by contra- 
diction. We select the first arc that succeeds after first failing and we exhibit a contradiction. Consider 
that arc (j,v ( - ) is this arc, which succeeded at stage g, but failed at least one earlier stage, and let / be the 
earliest such stage, f < g. It is straightforward to see that, in this case, ii < i% < ■ ■ ■ < if^\ < i g < if < 
if + i < • • • < ig-i, and arc (j,v; ) was tried and failed at all stages between / (inclusive) and g. 

As a thought experiment, let us stop the algorithm after step g. We have obtained g augmenting paths, 
thus a set P g , of g edge-disjoint paths, starting with arcs (s,vi), (s,V2), (s,Vf-i), (s,Vf), (j,v/+i), 

(s, v g ). Following the same though experiment, let us run the algorithm on digraph G', obtained from 
G, by deleting all arcs in V s except arc (s,V;) and those arcs that have been successful, before arc (s,vr) 
was first tried and failed. More formally, G' = (V,E'), where E' = E\ (V\V/), V' s = {0,vi), 0,v 2 ), 
(j,V/_i), (s,Vi )}. Obviously, \V^\ = f and digraph G' admits exactly / edge-disjoint paths, because (a) 
each of the remaining arcs in V' s can be the start of an edge-disjoint path in P g (which do not use any 
other arc in V s ), and (b) digraph G' cannot admit more than |V/| = / edge-disjoint paths. 

It is straightforward to see that, Algorithm [6] running on digraph G', will follow exactly the same 
steps as running on digraph G, up to the point when it first fails on arc (s, v, ) . At this point, the run on 
digraph G' stops, after finding / augmenting paths and constructing / edge-disjoint paths. 

Thus, the algorithm fails, because f < g, which contradicts its correctness. Therefore, the algorithm 
will never succeed on an arc that has already failed and never needs reconsidering again such arcs. 

This completes the second part of the proof. □ 

4 Structural and search digraphs in P systems 

In this section, we look at various way to reformulate the digraph edge- and node-disjoint path problems 
as a native P system problem. The P system we consider is "physically" based on a digraph, but this 
digraph is not necessarily the virtual search digraph G = (V,E), on which we intend to find edge- and 
node-disjoint paths. Given a simple P system IT = (0,K,8), where 8 is its structural digraph, we first 
identify cells as nodes of interest, V ~ K. However, after that, we see three fundamentally distinct sce- 
narios, which differ in the way how the forward and backward modes (i.e. backtrack and consolidation) 
of Algorithm [6] map to the residual arcs and finally to the structural arcs. 

1. We set E ~ 8. In this case, the forward mode follows the direction of parent-child arcs of 8, while 
the backward modes follow the reverse direction, from child to parent. 

2. We set E ~ {(v,u) \ (u,v) G 8}. In this case, the the backward modes follow the direction of 
parent-child arcs of 8, while the forward mode of the search follows the reverse direction, from 
child to parent. 

3. We set E ~ {(u,v),(v,u) \ (u,v) G 8}. In this case, the resulting search digraph is symmetric, 
and each of the arcs followed by the forward or backward modes of the search can be either a 
parent-child arc in the original 8 or its reverse. 

Cases (1) and (2) are simpler to develop. However, in this paper, we look for solutions in case (3), 
where all messages must be sent to all neighbors, parents and children together. Therefore, our rewriting 
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rules use the j8 and y G {repl} UK transfer operators (also indicated in Section]!]). Figure|7]illustrates 
a simple P module and these three scenarios. 




(a) (b) (c) 

Figure 7: Three virtual search digraphs for the same simple P module, (a) Same "physical" and search 
structure, (b) The search structure reverses the "physical" structure, (c) The search structure covers both 
the "physical" structure and its reverse. 

Note that, in any of the three cases, Algorithm [6] needs to be able to follow both the parent-child 
and the child-parent directions of P system structure. Therefore, the structural arcs must support duplex 
communication channels. 

After fixing the directions used by the virtual graph G, the next problem is to let the nodes identify 
their neighbors, i.e. discover the local network topology. 

5 Discovering cell neighbors 

In this phase, cells discover their own neighbors. Essentially, each cell sends its own ID to all its neigh- 
bors and records the IDs sent from its neighbors. This is a preliminary phase which is identical, for both 
edge-disjoint and node-disjoint versions. Table [8] illustrates the immediate neighborhoods, discovered 
at the end of this phase, for the P system of Figure [7] (a), with the virtual search structure shown in 
Figure [7] (c). 

Table 8: Neighbors table for the P system of Figure|7](a). An object nj indicates that cell a ; - is a neighbor 
of the current cell. 



Cell 


Neighbors 


Objects 


<?i 


{a 2 ,a 4 } 


{n 2 ,n 4 } 


o 2 


{ai, a 3 ,a 4 } 


{n u n 3 ,n A } 


°3 


{02,04,0 5 ,0 6 } 


{n 2 ,n 4 ,n 5 ,n 6 } 




{01,(72,(73,05} 


{ni,n 2 ,fi3,n 5 } 


o 5 


{a 3 ,a 4 ,a 6 } 


{n 3 ,n 4 ,n 6 } 


<?6 


{a 3 ,CT 5 } 


{m,n 5 } 



The set of objects used in this phase is {a, k,z,}U\Ji < j <n {g j, Uj,rij}. These objects have the following 
meanings: a indicates a cell reachable from a s ; k is the marker of the source cell; z is the marker of the 
target cell; tij indicates that a, is a neighbor of the current cell; gj, Uj indicate that Oj is the target cell; 
gj only appears in the source cell, while uj does not have this restriction. 
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Initially, the source cell a s has one copy of gj, representing the ID of the target cell Gj, and the other 
cells are empty. All cells start in state sq. Each reachable cell progresses through states sq,si,S2,st,,S4, 
according to the rules given below. In these generic rules (as elsewhere in this paper), we implicitly 
assume that (1) subscript i G {1,2, ...n} is customized for each cell to its cell ID; and (2) subscript j 
runs over all cell IDs (J G {1,2,.. .«}), effectively instantiating n versions of each generic rule where it 
appears. 

0. Rules for state 1- Rules for state si: 

1 so gi ->min so 1 si a -> min s 2 a (n/)$ rep i 

2 so gj ^min si ak 2. Rules for state s 2 : 

3 s„ «,- -+ min si az («,) w ! ^ fl ^ m . n 53 fl 

4 S() ^max Si 



5 so M; - -^ min si a (uj) u?1 



3. Rules for state st,: 

1 s 3 a -> min s 4 a 

2 *3 -> m ax S4 

The following example indicates how our generic rules are instantiated to take account the cell IDs, 
more specifically, how rules 0. 1 and 0.2 are instantiated in cell Oi : 

• SQ gl ->min SO 

• s gi -Wsi ak («i)^ repl 

• S gl ->min Si ak («2)^ repl 

• . . . , 

• s si ak (u n )$ 

The state transitions performed by cell a,-, i G { 1 , 2, . . . , n}, are briefly discussed below. 

• sq — > si : If Gi contains gj, Oi becomes the source cell. The source cell broadcasts object Uj to all 
its neighbors. After receiving an object uj, cell a,- becomes either the target cell, if i = j; or an 
intermediate cell, otherwise. Further, each cell relays one of received uj objects to all its neighbors. 

• si — > S2'. Cell Oi broadcasts to all its neighbors. Additionally, a,- accumulates nj objects from 
neighbors. 

• S2 — > S3: Cell Gj accumulates further n; objects from neighbors. 

• S3 — > S4. Cell Gi accumulates further nj objects from neighbors. Additionally, a, removes super- 
fluous ej objects. 

6 Simple P module rules for edge-disjoint paths algorithm 

First, we give a simple P module specification of the edge-disjoint paths algorithm presented in Section[3] 
We explicitly state our problem in terms of expected input and output. We need to compute a set of edge- 
disjoint paths of maximum cardinality between given source and target cells. 

Edge-Disjoint Paths Problem 

Input: A simple P module IT = (0,K,E, 8), where the source cell G s G K contains a token t t identifying 
the ID of the target cell G, G K. 

Output: If s ^ t, each cell a,- G K contains a set of objects P; = {pj | is a flow-arc} and a set of 
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objects C, = {cj | is a flow-arc} that represent a maximum set of edge-disjoint paths from o s to o t , 
where the following constraints hold: 

flow-arcs: c; G" Q, p t G" Pi, Cj G C, «4> p { G Py and c y - G Q => j £ 8(i) U 

source and target: P = and C, = 0. 

in flow = out flow: If i {s,?} then \Q\ = |P/|. 

only paths: With S(I) = \J ieI {j \ j G Q}, ^(I) = S(S(- ■ ■ 5(7) •••))= 0- 
Because of the network flow properties, we must also have \C S \ = \P t \, which also represents the maxi- 
mum number of edge-disjoint paths. 

This implementation has two phases: Phase I, which is the discovery phase described in Section [5] 
(using states so to 54), and Phase II, described below (which starts in state S4 and ends in state S13). Table[9] 
illustrates the expected algorithm output, for a simple P module with the cell structure corresponding to 
Figure [7] (a). For convenience, although these are deleted near the algorithm's end, we also list all the 
local neighborhood objects A^ ; - = {nj | j G 8(i) U (/)}, for i G {1,2, ... ,n}, which are determined in 
Phase I. 

Table 9: A representation of maximum edge-disjoint paths for simple P module of Figure|7](a). 



Cell \ Objects 


Ni 


Pi 


Q 









{C2,C 4 } 


o 2 


{721,723,714} 


{Pi} 


{^3} 


a 3 


{?22, 72 4 , 72 5 , 72 6 } 


{P2,Pa} 


{C5,C 6 } 


a 4 


{721,722,723,725} 


{pi} 


{C3} 


05 


{723, 72 4 , 72 6 } 




{C6} 


(76 


{m,n 5 } 


{^3,^5} 






In addition to the set of objects used in Phase I, this phase uses the following set of objects: {bj, Cj, 
dj, ej, fj, hj, nij, pj, qj, rj, tj, Xj, y ; } U {v, w}. In cell a,-, these objects have the following meanings: 

• bj indicates a pushback, received from a,'s flow-successor oy, 

• ej records the sender of a pushback, if a,- is not-yet-visited; 

• hj records the sender of a pushback, if a,- has already been visited; 

• Cj indicates that Oj is a,'s flow-successor; 

• pj indicates that Oj is (7/'s flow -predecessor; 

• rj records a pushback sent by a,- to its flow-predecessor a,-; 

• tj records a backtrack request, after a failed pushback to (Jj\ 

• dj indicates that Oj is (7,'s search-successor; 

• qj indicates that Oj is a,'s search-predecessor; 

• fj indicates an attempted search extension received from Oy; 

• itij records that Oj was unsuccessfully tried; 

• Xj indicates that Oj rejected a flow-extension or flow-pushback attempt; 
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• Xj indicates that Oj accepted a flow-extension or flow-pushback attempt; 

• v requests that a, resets the record of all tried and visited neighbors; 

• w requests that a, remains idle for one step. 

Initially, the source cell o s has one copy of gj, representing the ID of the target cell Oy, and the 
other cells are empty. All cells start in state 5o. According to the rules of Phase I, each reachable cell 
progresses to state 54, which is the start of Phase II, whose generic rules are given below. In these rules (as 
in Phase I), we implicitly assume that (1) subscript / G {1,2, ...n} is customized for each cell to its cell 
ID; and (2) subscripts j,k run over all cell IDs, j,k £ {1,2, . . .n}, and j 7^ k. To apply the optimization 
proposed in Section 3.4 replace rule 5.3 "55 djXj — hnin ^5 am" by "55 djXj - 



4. Rules for a cell a, in state 54: 

1 s 4 k -w s 5 k 

2 54 Z — ?"min ^6 Z 

3 S4 a — hnin ^7 a 

5. Rules for a cell a, in state 55: 

1 s 5 anj dj (fi)$. 

2 s 5 d jyj ^ fflin s l2 acjww (v)^ repl 

3 s 5 djXj -^ min s 5 anij 

4 s 5 bj -W5 5 (xi)$. 

5 s 5 fj -> min s 5 (xi)$j 

6 s 5 ak -w ^13 aaww (a)^ repl 

6. Rules for a cell a, in state s$: 

1 s 6 njfj ^ min s 6 pj (yi)$. 

2 s 6 v -> min s i2 ww (v)^ repl 

3 5 6 aaz -^min ^13 <www (a)^ repl 

7. Rules for a cell a, in state 57: 

1 $7 V -> mln ^12 WW (v)^ repl 

2 57 -> min 513 aaww (a)^ repl 

3 57 njfj 

4 57 cyfcy -> min 5 8 ey 

5 5 7 /zy -W^ll Cj (Xi)$. 

6 5 7 J9y^ -> m in^lO Pjqk 

7 5 7 <?y -W5 7 my (x,)^. 

8 ^7 fj ->min ^7 (jCf)^ 

8. Rules for a cell a,- in state s&: 

1 5 8 anj -» min 59 a<iy (/ ; )^ 

2 5 8 a -> min 5 K ) a 



n 55 a . 

9. Rules for a cell a, in state 59: 

1 59 djyjet -w 57 Cjin k (y,-)^ 

2 59 dyyytf* -> min s 7 Cjp k (y,-)^ 

3 59 dyXy ~W 5 8 my 

4 5 9 Cy6y -> min 5 9 Cj (Xi)^. 

5 59 ray/y -^ m in s 9 rnj [xi\. 

10. Rules for a cell a, in state 510: 

1 510 apj -> min 5n arj (fe/)^ 

2 510 aey 5 7 OCy (x,-)^. 

3 ^10 y ^min 

57 anij (xi)$. 

11. Rules for a cell a, in state s\\ : 

1 5n 0'yj e * "^min S7 mym* (y,)^ 

2 5n ^jft ~W S7 WyPA (ji)^ 

3 511 O'.V/ > :r .ir: 5||) // 

4 5n Cyfey -WS7 /iy 

5 5n «y/; -W *11 rnj {Xi\, 

12. Rules for a cell a,- in state 512: 

1 S U W -W 512 



2 5i2 V — 

3 5i2 my 

4 512 tj - 

5 5i2 £ — 

6 5i2 Z - 

7 512 « 



max ^ 12 



■>min *12 
^min ^5 ^ 

'min ^7 a 

13. Rules for a cell a, in state 513: 

1 *13 W-^ min 5i3 

2 513 a — ^ max *0 

3 513 tj -Wi *0 

4 513 «/-WJo 

5 5 i3 nij -^min^O 
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The following paragraphs describe details of several critical steps of our edge-disjoint algorithm, 
such as forward and consolidation modes of intermediate cells. 

• The rules in state s% cover the forward mode attempt of an intermediate cell a,, via a non-flow arc. 

If Oj has a not-yet-tried neighbor Oh $ PiUQ, the search continues with <j/,; otherwise the search 
backtracks. 

• The rules in state sg cover the consolidation mode for an intermediate cell a,, who has succeeded 
a forward extension on a non-flow arc. 

During the consolidation process, the behavior of a, depends on whether a, was reached by a flow 
arc or non-flow arc: 

o If Oi was reached by flow arc in a reverse direction, a, replaces its flow-successor Oj 
with Oh, where ofc is its search-predecessor. 

o If Oi was reached by a non-flow arc (k,i), a, sets Cfy as a flow-predecessor and Oh as a flow- 
successor, where Oh and Oj are a,'s search-predecessor and search-successor, respectively. 

• The rules in state sio cover the forward mode attempt of an intermediate cell a,, via a pushback. 

If Oi has a not- yet-tried flow-predecessor O/,, the search continues with Oh (i.e. use flow arc (h,i), 
in a reverse direction); otherwise the search backtracks. 

• The rules in state su cover the consolidation mode for an intermediate cell a,, who has succeeded 
a forward extension via a pushback to Oh (i.e. on a flow-arc (h, i), in a reverse direction). 

During the consolidation process, the behavior of a, depends on whether a, was reached by a flow 
arc or non-flow arc: 

o If Oi was reached by flow arc (i, j), in a reverse direction, a, removes its flow -predecessor Oj, 
and flow-successor Oj. 

o If Oi was reached by a non-flow arc (k,i), Oi replaces its flow -predecessor G), with Ok, where 
Ok is (Tj's search-predecessor. 

Theorem 9. For a simple P module with n cells and m=\8\ edges, the algorithm in this section runs in 
0{mn) steps. 

7 Simple P module rules for node-disjoint paths algorithm 

First, we give a simple P module specification of the node-disjoint paths algorithm presented in Section|3] 
We explicitly state our problem in terms of expected input and output. We need to compute a set of node- 
disjoint paths of maximum cardinality between given source and target cells. 

Node-Disjoint Paths Problem 

Input: A simple P module n = (0,K,E, 8), where the source cell o s € K contains a token t t identifying 
the ID of the target cell o, € K. 

Output: If s ^ t, each cell a,- G K contains a set of objects P, = {pj \ (j,i) is a flow-arc} and a set of 
objects Cj = {cj | (i,j) is a flow-arc} that represent a maximum set of node-disjoint paths from o s to o t 
where the following constraints hold: 
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flow-arcs: q Q, p t Pj, cj G Q 44> G P; and c ; - G C, =^ 7 G 5(i) U <5 _1 (/). 
source and target: P s = and C, = 0. 
node-disjoint: If i G" {s,t} then |C,-| = |P,-| < 1. 

onlypaths: WithS(/) = { *. ^^^f = ° (0 = W • •■•))= 

Because of the network flow properties, we must also have \C S \ = \P t \, which also represents the maxi- 
mum number of node-disjoint paths. Notice the constraints to require only paths has been simplified in 
that the successor S(i) of non-source cell a, is a single cell instead of a set of cells that was needed for 
the general edge-disjoint problem. 



Table 10 illustrates the expected algorithm output, for a simple P module with the cell structure 
corresponding to Figure [7J (a). For convenience, although these are deleted near the algorithm's end, we 
also list all the local neighborhood objects Ni = {nj \ j G 8(i) U (/)}, for i G {1,2, .. . ,«}, which are 
determined in Phase I. 



Table 10: A representation of maximum node-disjoint paths for simple P module of Figure[7](a). 



Cell \ Objects 


N 


Pi 


Q 




{n 2 ,n 3 } 





{c 2 ,c 4 } 


O2 


{721,723,714} 


{Pi} 




O3 


{?22,724,72 5 ,72 6 } 


{P2} 


{<*} 


O4 


{72i, 72 2 , 723 , 72 5 } 


{Pi} 




CT5 


{7I3,72 4 ,72 6 } 






<?6 


{m,n 5 } 


{P3,Ps} 






The rules of this node-disjoint path algorithm are exactly the rules of the edge-disjoint path algorithm 
described in Section[6| where the rules for state s 7 are replaced by the following group of rules. 

7. Rules for a cell a, in state s 7 : 

1 s 7 v ^ min s n ww (v)^ repl 

2 s 7 aa -W S13 aaww 0)^ repl 

3 s 7 njfjp k 

hnin ^10 9 jPk 

4 s 7 njfj 

~ ^min s 8 Q j 

5 s 7 cjbj -Ws 8 ej 

6 sjhj -w ^8 ej 

1 S 7 ry-)- min S\\ Tj 

o s 7 qj 7- m i n s 7 rrij (xi)$j 
9 s 7 fj -> mln s 7 (xi)$. 

The new state s 7 rules implement our proposed non-standard technique described in Section 3.2 for 
enforcing node capacities to one, without node-splitting. 

The running time of our node-disjoint paths algorithm runs in polynomial number of steps, since the 
algorithm is direct implementation of Ford-Fulkerson's network flow algorithm. 

Theorem 10. For a simple P module with n cells and m=\8\ edges, the algorithm in this section runs 
in 0(mn) steps. 
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8 Conclusion and Open Problems 

Using the newly introduced simple P modules framework, we have presented native P system versions 
of the edge- and node-disjoint paths problems. We have started from standard network flow ideas, with 
additional constraints required by our model, e.g., cells that start without any knowledge about the local 
and global structure. Our P algorithms use a depth-first search technique and iteratively build routing 
tables, until they find the maximum number of disjoint paths; Our P algorithms run in polynomial time, 
comparable to the standard versions of the Ford-Fulkerson algorithms. We have proved and used a 
speedup optimization, which probably was not previously known. For node-disjoint paths, we proposed 
an alternate set of search rules, which can be used for other synchronous network models, where, as in 
P systems, the standard node-splitting technique is not applicable. 

All of our previous P algorithms assumed that the structural relation 8 (of a simple P module) sup- 
ports duplex communication channels between adjacent P system cells. Substantial modifications are 
needed when we consider the simplex communication case. It is not just a simple matter of changing the 
rules of the systems to only following out-neighbors when we are finding paths from the source to the 
target — we have explicitly utilized the ability to "push back" flow on a flow-arc, by sending objects to 
their flow-predecessors, when hunting for augmenting paths. Thus, some new ideas are needed before 
we can compute disjoint paths when the structural arcs allow only simplex communication. 

We also want to know if we can we solve the problem of finding disjoint paths between k pairs of 
(s\,t\) . . . (sic^ic), that is comparable in performance to the 0(n 3 ) algorithm of [16]. 

We are interested to know whether a breadth-first search (BFS) approach would be more beneficial 
than a depth-first search. By using BFS we could potentially exploit more of the parallel nature of 
P systems. 

By combining this paper's results with our previous P solutions for the Byzantine problem [3j|2l, 
we have now solved one of our original goals, i.e. to solve, where possible, the Byzantine problem for 
P systems based on general digraphs, not necessarily complete. This more general problem can be solved 
in two phases: (1) determine all node-disjoint paths in a digraph, assuming that, in this phase, there are 
no faults; and then (2) solve the consensus problem, even if, in this phase, some nodes fail in arbitrary, 
Byzantine ways. An interesting problem arises, which, apparently, hasn't been considered yet. What can 
we do if the Byzantine nodes already behave in a Byzantine manner in phase (1), while we attempt to 
build the node-disjoint paths? Can we still determine all or a sufficient number of node-disjoint paths, in 
the presence of Byzantine faults? 

Some of previous experiences [UQI3] have suggested that P systems need to be extended with support 
for mobile arcs, to incrementally build direct communication channels between originally distant cells 
and we have offered a preliminary solution. The current experience suggests that this support should 
be extended to enable straightforward creation of virtual search digraphs on top of existing physical 
digraphs. 

We consider that this continued experience will provide good feedback on the usability of P sys- 
tems as a formal model for parallel and distributed computing and suggest a range of extensions and 
improvements, both at the conceptual level and for practical implementations. 
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Table 11: Edge-disjoint paths solution traces (steps 0,1,... ,30) of the simple P module shown in Fig- 
ure [7] (a), where 0"i is the source cell and G(, is the target cell. 



Step\Cell 




o 2 


<7 3 


04 


0"5 


0-6 





J g6 


so 


■SO 


so 


•SO 


■SO 


1 


si ak 


so U 6 


i'O 


sq u 6 


4"0 


■SO 


2 


S2 alaii 


s\ an\U(, 


J Kg 


s\ an\u^ 


JO 1(6 


■SO 


3 


J3 aknin\iR 


2 

S2 an\n 4 u. 


2 

j'l 0/72/74/75 


3 

J2 an\n2U^ 


Jl 0/74/75 


J Mg 


4 


S4 akn^rii, 


2 

J3 an\nT,n 4 u^ 


3 

J2 an2n\ne ! u^ 




2 

J2 an3«4Mg 


Jl 0/73/752 


5 


J5 akti2ri4 




3 

J3 av\2n\n^n^u b 






2 

J3 0/73/74775/7:: 


J2 0/73/752 


J4 ani«3n 4 


s 4 ani^nyis 






















6 


j'5 d\kt\2 


J7 a/71773774 


J4 0772/74/75/75 




si a/imn2»3«5 




J4 0/73/74/75 


J3 0/73/752 




7 

i 


55 d 4 kii2 


J7 a/71/73774 


Si an2n\n^n^ 


jg an 2 n 3 n 5 q ] 


J7 £0/73/74/75 


J 4 0/73/752 




Q 
O 


J5 d 4 kii2 


s-j (in\Yi3n\ 


si af 4 n 2 n 4 n 5 n 6 


S9 arf3«2' 7 5?l 


J7 a/73/74/75 


J5 0/73/752 


Q 


j'5 d 4 kii2 


J7 CLYl\H r i } Yl 4 


s s an 2 n 5 n 6 q 4 


J9 ad3«2 f! 5?l 


J'7 a/73/74775 


J5 0/73/752 


1U 


J 5 d 4 kii2 


J7 a/3 77 1/73 7/4 


J9 ad2n^n(,q 4 


S9 arf3«2' 7 5?l 


J7 a/73/74/75 


J5 0/73/752 


i i 


55 d\kvi2 


Jg Gfl\fl^C/3 


s 9 ad 2 n 5 n 6 q 4 


S9 003772/7501 


J'7 a/73/74775 


J5 0/13/752 


12 


Js ad fokno 


J9 ad\n 4 q-$ 


J9 adjn^n^q^ 


s<) adj,n2n^q\ 


J7 a/73/74/75 


J5 0/73/752 


13 


55 d\ki\2 


Sg ad\n 4 q-x,x\ 


j 9 ad 2 n 5 n 6 q 4 


s<) adj,n2ii^q\ 


J7 fl/73/74/75 


J5 0/73/752 


14 


J5 04/072 


Jg £JOTl«4^3 


J9 adjn^n^q^ 


s<) adj,n2n^q\ 


J7 a/73/74/75 


J5 0/73/752 


15 


j'5 d\kn.2 


J9 ad/\m\q3 


j'9 002/75/7504 


sg ad3f 2 n 2 n 5 qi 


J7 0/737/4/75 


J6 0/13/752 


16 


J5 d\kii2 


sg 004/77103x4 


J9 002/7577504 


j'9 003/772/7501 


J7 a/73/74/75 


J5 0/73/752 


17 


J5 d\kri2 


jg 0/771/77403 


sg ad2n^n(,q 4 


sg ad3in2ii^qi 


J7 a/73/74/75 


J5 0/73/752 


18 


j'5 d\ki\2 


S\o am\m 4 qj, 


J9 002/75/7504 


sg ad3rri2n^qi 


J7 a/73/74/75 


J5 0/73/752 


19 


J5 d\kn.2 


sj am\m-},m\ 


sg ad2n$n(,q 4 X2 


sg 003/772/7501 


J7 a/73/74/75 


J5 0/73/752 


20 


j'5 d\kn.2 


Si amim^m^ 


J8 07772/75/7504 


sg ad3in2n^qi 


J7 a/73/74/75 


J 6 0/73/752 


21 


J5 d\kr\2 


s-j am\m-},m\ 


J9 005/772/7504 


sg 003/772/7501 


J7 0/73/74/75 


J 6 0/3/73/752 


22 


j'5 d 4 kri2 


si am\mj,m 4 


J9 ad 6 m2n 5 q 4 y 6 


sg ad3rri2n^qi 


J7 0/737/4/75 


J 6 0/75/332 




23 


J5 d\kn2 


S7 aniim^m^ 


si ac 6 m 2 n 5 p 4 




sg adii^nsqiyi 


J7 a/73/74/75 


J6 0/75/732 


24 


J 5 rf 4 te 2 )'4 


S7 amim^m^ 


Si ac(,m2ri5p 4 


Si aC2ni2n^pi 




J7 a/73/74/75 


J 6 0/75/732 


25 


Jl2 ac4fcn2>v 2 


si 0/771/773/774 v 


si ac(,m2ri5p 4 


si acT,m2ns,p\v 


J7 0/73/74/75 


J 6 0/75/732 


26 


j"12 ac 4 £/72V 2 w 


2 

j 12 am\m-itna,wr 


Si aC(,m2ri5p 4 v 


I 


2 

J 12 af3m2«5Pi vw 


J7 0/73/74/75V 


J 6 a/7 5 p 3 2 


27 


J12 ac 4 kri2 


s\2 an\n-},n 4 vw 


j'12 ac 6 /772«5P4V 2 H' 2 


j'12 acj,n 2 n 5 piv 




2 

J12 0/73/74/75WV 


j 6 a/7 5 p 3 v 2 2 


28 


J5 aci\kn2 


s\2 an\ri3n 4 


J12 flC6«2«5P4i'H' 


s i2 ac 3 n 2 n 5 pi 


J12 a/73/74/75 VW 


J12 a/75P3V>V 2 2 


29 


J5 c\d2k 


Si afin\n%n 4 


J12 aC6«2«5P4 


si ac3H2ns > p\ 


Jl2 a/73/74/75 


J12 a/7 5 // 3 W2 


30 


J5 c\d2k 


s% 0/73/7401 


J7 ac^ii2n^p 4 


si ac^U2n 5 p\ 


J7 a/73/74/75 


J12 a/7 5 p 3 2 
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Table 12: Edge-disjoint paths solution traces (steps 31,32, ... ,61) of the simple P module shown in 
Figure [7] (a), where di is the source cell and 0$ is the target cell. 



Step\Cell 


<5\ 


02 


CT 3 


04 






31 


s 5 c 4 d 2 k 


.59 adi l ri4q\ 


«7 ac 6 f 2 n 2 n 5 p4 


Si acj,n 2 n^pi 


S7 a«3n4«6 


56 an^p^z 


32 


S5 c 4 d 2 k 


sg ad^U4q\ 


s 8 ac 6 n 5 p 4 q 2 


Si acj,n 2 n^pi 


S7 fln3n4«5 


Sg flW 5 p3Z 


33 


s 5 c 4 d 2 k 


sg ad^ri4q\ 


Sg ac 6 d 5 p4q 2 


s 7 aci,n 2 n 5 pi 


s 7 af 3 n 3 n 4 n 6 


56 anspj,z 


34 


S5 c 4 d 2 k 


sg ad^H4q\ 


Sg ac 6 d 5 p 4 q 2 


Si ac 3 n 2 n 5 p l 


sg an 4 n 6 q 3 


sg arij-piz 


35 


s 5 c 4 d 2 k 


sg ady%4q\ 


Sg ac 6 d 5 p 4 q 2 


si ac 3 f 5 n 2 n 5 pi 


sg arf4«g^3 


56 an^p^z 


36 


S5 c 4 d 2 k 


sg ad^ri4qi 


Sg ac 6 d 5 p4q 2 


sg ac 3 n 2 Pl "75 


S9 ad4ri( J q-i, 


56 ansp^z 


J l 


Sg c 4 d 2 k 


Sg Cld^f4Vl4q\ 


sg ac(,dsP4q 2 


sg acT,d 2 piq$ 


sg cid4ri( J qT l 


56 CIII5P3Z 


38 


55 c 4 d 2 k 


sg cid^m4q\ 


sg ac&ds V4Q 2 


sq ac\d~> v\ (7sxo 


SQ Culdtlf,QT. 


Sf. an^ i?rtz 


39 


55 c 4 d 2 k 


sg ad^m.4q\ 


Sg ac 6 d 5 p4q 2 


sg acyn 2 p\q 5 


S9 ad4ti( > q3 


56 an 5 pi,z 


40 


55 c\d 2 k 


sg ad-^m4q\ 


sg ac 6 d 5 p 4 q 2 


Sio ac 3 m 2 pi?5 


sg ad4fi( I qT > 


s 6 a«5P3Z 


41 


55 bi\c\d 2 k 


sg ad^m.4q\ 


sg ac 6 d 5 p 4 q 2 


su acj,m 2 q 5 r\ 


sg ad4n(,q-i, 


56 a«5P3Z 


42 


S5 ce,d 2 k 


sg ad^m,4q\ 


sg ac 6 d 5 p 4 q 2 


su acT,m 2 q$r\x\ 


sg ad4ri(,qj, 


s 6 a«5P3Z 


43 


55 c 4 d 2 k 


sg ad^m.4q\ 


sg ac 6 d 5 p 4 q 2 


Sio acim 2 q 5 ti 


sg ad4ti( J q3 


56 ««5P3Z 


44 


55 c 4 d 2 k 


sg ad^m,4q\ 


sg ac 6 d 5 p 4 q 2 


S7 ac},m 2 m^t\ 


sg ad4ri(,qj,x 4 


s 6 a«5P3Z 


45 


55 c\d 2 k 


sg ady,ni4qi 


Sg ac^d 5 p4q 2 


si acjmzmsH 


sg flm4«5^3 


56 ««5P3Z 


46 


55 c\d 2 k 


sg ad^m4q\ 


sg ac 6 d 5 p 4 q 2 


si acT,m 2 ms,ti 


S9 ad(,ni4q2 


s 6 af 5 n 5 p 3 z 


47 


s 5 cs,d 2 k 


sg ad^ni4qi 


Sg ac 6 d 5 p 4 q 2 


si ac^m 2 m^t\ 


Sg ad 6 m4q 3 y b 


*6 «P3P5Z 


48 


55 CA,d 2 k 


sg ady,ni4qi 


sg ac 6 d 5 p4q 2 y 5 


si ac^in 2 m$t\ 


si ac(,m4P2 


S6 "P3P5Z 


49 


s 5 Ci,d 2 k 


sg ady l ni4qiy2 


s 7 ac 5 c 6 p 2 p 4 




Si ac^m 2 m^t\ 


Si aC(,tri4P3 


56 apip^z 


JU 


*5 c \d 2 ky 2 


Sj CLC^}714P\ 


s^ ac 5 c 6 p 2 p 4 


si cict ) Yn 2 vtiz j t\ 


Si aC(,m4P2 


56 apipsz 


J 1 


s\ 2 


57 WI4 P\V 


s 7 ac 5 c 6 p 2 p4 


Cr~ ! 1 1 ' 1 111 ~. Ill ^ / , 1 1 

57 ClL^llL 2 rtl^l\V 


Si cic , £ ) ffi4P3 


s'6 apT,p$z 


52 


s\ 2 ac 2 C4kv 2 w 


2 

s\2 aCT,m4p\vw 


S7 ac 5 c 6 p 2 p 4 v 


2 


S\ 2 ac^m 2 m^t\vw 2 


Si ac b m4p^v 


5 6 ap3P5Z 


53 


S\ 2 ac 2 c\k 


s\ 2 ac2ti4pivw 


*12 ac 5 c 6 p 2 p 4 v 2 w 2 


s\ 2 acj,n 2 n^p\v 2 w 


s\ 2 ac^ni4p^vw 2 


s 6 ap 3 p 5 v 2 z 


54 


55 ac 2 Ci\k 


s u ac 3 n 4 pi 


■s 12 ac 5 c 6 p 2 p 4 vw 


S12 ac 3 n 2 «5Pl 


S12 ac6/34/?3Vw 


S12 ap3P5VW 2 z 


55 


S13 a 2 c 2 C4\v 2 


57 a 2 cj,ri4p\ 


*12 ac 5 C(,p2P4 


2 

S7 a C3«2 ra 5Pl 


S12 ac 6 n 4 p 3 


S12 ap 2 p 5 wz 


56 


S13 a A c 2 c\w 


S13 a i c^n 4 piw 2 


s 7 a 3 c 5 c 6 p 2 p4 


S13 a 3 C3W2fl5PiW 2 


si a 2 c 6 ii4pi, 


H2 aPiPiZ 


57 


4 

S13 a c 2 C4 


4 

S13 a C3M4P1W 


«13 a 4 c 5 C(,p 2 p4W 2 


S13 a^C3«2 n 5/'l> 1 ' 


3 2 
S13 a cg«4P3W 


s 6 a 3 pip 5 z 


58 


so c 2 c 4 


513 a 4 C3«4Pi 


S13 a 5 c 5 c 6 p 2 p 4 w 


S13 a 5 cin 2 n 5 p\ 


4 

S'13 a C6ii4P3W 


S13 a 3 P3P5W 2 


59 


s c 2 c 4 


SQ CiPl 


513 a 5 c 5 c 6 p 2 p 4 


so C3P1 


513 a 4 C 6 «4P3 


S13 « 3 P3P5>v 


60 


so c 2 C4 


s cipi 


■so c 5 c 6 p 2 p4 


so C3P1 


so C6P3 


513 a 3 P3P5 


















61 


so c 2 C4 


S c 3 p! 


SO CsC 6 p 2 p4 




so C3.Pl 


SQ C6P3 


■SO P3P5 
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Table 13: Node-disjoint paths solution traces (steps 0,1,... ,29) of the simple P module shown in Fig- 
ure [7] (a), where 0\ is the source cell and o& is the target cell. 



Step\Cell 


°1 


a 2 


CT3 


<T 4 


C7 5 


0-6 





JQ gg 


JQ 


•so 


SO 


JQ 


S() 


1 


Jl a£ 


JO 


Jo 


SQ u 6 


J() 


■SO 


2 




Ji a/21 


J Wg 


S{ an\U(, 


JO Z-/g 


so 


3 


^ j v Z M- n 


2 

z 1 ^+ ll 


2 

Jl aii2«4Mg 


3 

J2 an\ii2W^ 


jj aii4U(y 


J 


A 




2 

J3 (lll\tlT l H^ll^ 


3 

J2 a/?2' ! 4' ! 5"g 


3 

J3 a;3i«2"3' I 5 i 'g 


2 

S2 £1/23 '^4^g 


Jl a«3«5Z 






















5 




J4 flni«3M4 


J3 flft2' I 4' I 5''6 t '(i 


J4 aMin2«3«5 




2 

J3 an^n4n^uf^ 


J2 an^n^z 
























6 


^5 d<\kri2 


J7 anxn^n^ 


J4 fl«2«4"5' 7 6 




J7 a/i«in2' 3 3«5 




S4 an^n4ii(, 


J3 an^n^z 




7 




sj an\n^n\ 


J7 a«2' ! 4' J 5"6 


Jg flW2«3«5?l 


J7 an3«4«6 


S4 an^n^z 




8 




J7 ani«3n4 


J 7 a/ 4 n2«4«5«6 


J9 adi^n^qi 


J7 an^^nf, 


jg a«3)i 5 z 


9 


^5 d^kn2 


J7 an\n^n\ 


Jg an2' ! 5' ! 6 l ?4 


J9 arf3«2"5?l 


J7 an3;34«6 


jg ani,n 5 z 


10 


s$ d^kni 


J7 afsninyity 


J9 ad2n 5 n 6 q4 


J9 adi^n^qi 


J7 anj,n4n( > 


Jg a«3)lgz 


11 


^5 d^kri2 


jg a«i«4(73 


sg ad2n^n(,q4 


J9 adj,n2n^qi 


J7 an3n4«6 


j 6 a« 3 M 5 z 


12 


s 5 d 4 f 2 kn 2 


jg ad\nt\qs 


sg ad2n^ii(,q4 


J9 adj,n2n^q\ 


J7 an^^rib 


jg anyi^z 


13 


s$ d 4 kri2 


J9 arfin4</3Xi 


J9 ad2n^n(,q4 


J9 adj,ri2n^q\ 


s-j an^n4ri(, 


jg an^n^z 


in- 


5*5 d^kn.2 


Jg if /77 1 724^/3 


J9 ad2ns l n( l q l \ 


J9 ad^n2n=,q\ 


Si ci/73'^4'^6 


jg a«3ii 5 z 


15 


ddkno 


J9 ad^m\q^ 


J9 ad2n^ri( l q l \ 


j 9 ad3f 2 n 2 n 5 qi 


J7 a/73/74/76 


jg a«3(? 5 z 


16 


5< ddkno 


Sg ad^niiq-^x^ 


J9 ad2"5' J 694 


J9 ady,rri2n^qi 


J7 (1113114116 


Jg fl«3)l 5 z 


17 


ddkno 


jg am\m\q^ 


J9 ad2n$n( l q l \ 


J9 adT,rri2n^q\ 


si (111^411^ 


jg a«3« 5 z 


18 


■Ss ddkno 


Jin am\m±Qx 


J9 ad2ns l n( l q l \ 


J9 adT,ni2n^qi 


si (112311411^ 


Jg fl«3)1 5 z 


19 


Js ddkni 

J H Z 


J7 aniim^m^ 


J9 ad2iis l n( l q l \X2 


J9 adT,rri2n^qi 


Si (11x311411^ 


jg an^n^z 


20 


J'5 d^kn^ 


J7 am\mj,m^ 


Jg am2 n 5 n 6<?4 


J9 ady,rri2n^qi 


si aniens 


jg an.T,n 5 z 


21 


J5 d^ktl2 


J7 am,\m^m\ 


J9 ad(,m2ii^q\ 


sg adx,ni2n^q[ 


Si an3;34«6 


jg af 3 n 3 n 5 z 


22 


j'5 d\ki\2 


J7 am \ 17131114 


s 9 ad 6 m 2 n 5 q 4 y 6 


sg adT,rri2n^qi 


J7 anj,n4n(, 


s 6 an 5 p 3 z 




23 


j'5 d\ki\2 


J7 am\mj,m^ 


s^ ac^rri2n 5 p4 




sg adT,m2n^q\y3 


Si anj,n4ti( > 


jg an 5 p 3 z 


24 


j 5 d 4 kn 2 y4 


J7 am,\m^m\ 


J7 ac(,ni2n^p4 


j 7 ac 3 m 2 n 5 p 1 




Si an3«4n6 


j 6 a« 5 p 3 z 


25 


s\2 acA,kniw 2 ' 


J7 amim^m^v 


s^ ac^ni2nsP4 


Sj ac^ni2n^piv 


J7 an3«4n6 


Jg a«5P3Z 


26 


s\i ac4kn2V 2 w 


S\2 amxm^m^vw 2 


J7 acgm2«5P4V 


Z 


Jl2 ac3m2n5Pivw 2 


J7 an3«4«6V 


Jg a«5P3Z 


27 


S\2 ac\kri2 


s\2 an\nj,riA.vw 


2 2 

Jl2 aC6m2«5P4V w 


j I2 ac 3 n2«5Piv 


2 

w 


2 

j]2 aH2n4ii(,vw 


jg an^p2v 2 z 


28 


J5 ac\kn2 


s\2 an\riT,n.i\ 


Jl2 flC6n2"5P4 vvv 


Jl2 ac 3 n 2 n 5 pi 


Jl2 a«3fi4n6VW 


S\2 ansp3vw 2 z 


29 


j'5 c\d2k 


J7 af\n\n^,ni\ 


Jl2 ac 6 n 2 n 5 p 4 


s-j acj,n2ns,p\ 


Jl2 fl«3«4«6 


Jl2 an5P3Wz 
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Table 14: node-disjoint paths solution traces (steps 30,31, ... ,59) of the simple P module shown in 
Figure [7] (a), where di is the source cell and Of, is the target cell. 



Step\Cell 




0-2 


(*3 


04 






30 




Sg an^iiA.qi 


sj ac 6 n 2 n 5 p4 


sj ac^nin^pi 


*7 aniens 


S12 a«5P3Z 


31 


S5 ci\djk 


sg ad 3 n 4 q\ 


sj ac 6 f2n 2 n 5 p4 


sj ac$n2ti^pi 


sj an3«4n6 


S 6 a«5f 3Z 


32 


55 CA,djk 


sg ad^ntyqx 


•S10 ac 6 n 5 p 4 q 2 


sj ac3ri2n^pi 


i'7 an-$n4n(, 


S 6 an 5 p3Z 


33 


S5 Ct\djk 


sg ad^n^qi 


ill ac 6 n 5 q 2 r4 


sj ab-i,c^n2n%p\ 


s-j an-t,Yi4n^ 


sg a«5P3z 


34 


s 5 c 4 d 2 k 


sg ad^ntyqx 


Su ac 6 n 5 q 2 r4 


s% ae 3 n%ns,p\ 


sj an^n4n(, 


S 6 a«5/?3Z 


35 


S5 Ci\djk 


Sg ad-},fA,n\q\ 


ill ac 6 n 5 q 2 r4 


sg ad2ej,n^p\ 


sj an^n4ii(, 


S 6 a«5P3Z 


36 


S5 C4Cfefc 


sg ad^m^qi 


Su ac 6 n 5 q 2 r4 


sg ad2e^,n^p\X2 


sj ari2n4ti(, 


S 6 a«5/?3Z 


37 


£5 c^d^k 


sg dd 3 n\4q\ 


•Ml LtL 6 /i 5 l jz' 4 


•>& wc 3" t z' t 5i / l 


Sj ClH^H4flfo 




38 


S5 C4fi?2^ 


sg ad 3 m.4q\ 


ill ac 6 n 5 q 2 r4 


sg ad^ej,fri2P\ 


Sj a] 4n3n4«5 


s 6 an 5 p 3 z 


39 


S5 C\d2k 


sg ad^mi\q\ 


ill ac 6 n 5 q 2 r4 


sg ad^e^m2P\ 


Sg ailing 


s b an 5 p 3 z 


40 


55 CA,djk 


sg ad^m^qi 


su ac 6 f 5 n 5 q 2 H 


sg ad^eo,m2P\ 


sg ad-$ii( l q4 


S 6 a«5/?3Z 


41 


S5 Ct\djk 


sg ad^m^qi 


su ac( s mz,q2r4 


sg ads > e r i,m2P\ 


Sg fld3«g^4X3 


S6 a«5P3Z 


42 


S5 C4fi?2^ 


sg ad^m^qi 


su ac 6 m 5 q 2 r4 


sg ad$ej,m2P\ 


Sg amj l n( i q4 


S 6 a«5/?3Z 


43 


55 Ct\djk 


sg ad^m^qi 


su ac b m i q 1 r4 


sg ad^e^m2P\ 


sg ad(yinj,q4 


■S6 af 5 n 5 p3Z 


44 


S 5 c 4 d 2 £ 


sg ad^m^qi 


s\\ ac 6 m 5 q 2 r4 


sg adsej,m2P\ 


sg ad b m 3 q4y(, 


*6 "P3P5Z 


45 


^5 c^d^k 


sg cid^vyi^qi 




vn fldzPiffln r>\ v< 


Sj QCfcm2P4 


V£ /7 771 77<: 7 


A £. 

40 


J 7. 

S5 c 4 d 2 K 


sg ad^niiyqi 


ill ac 6 m 5 q 2 r4y4 


i 7 ac 5 m 2 m3Pi 


Sj aC(,m^p4 


■S6 "P3P5Z 


47 


55 Ct\djk 


sg adynnqiyj, 


sj ac(,m4m$p2 


57 ac5iii2>n3p\ 


Sj ac(,m^p4 


■S6 «P3P5Z 


A Q 


S5 C4d2ty2 


sj acpn^pi 


sj ac( J m4m^p2 


Sj ac$m2m 3 p\ 


sj ac( J m^p4 


s 6 ap 3 p 5 z 


49 


s 12 ac^CAkw 1 


sj ac^m^p jv 


sj ac( 1 m4ins t p2 


Sj ac^m2in^p\v 


sj ac b m 3 p 4 


■S6 «P3P5Z 


50 


s\2 ac2C\k\^w 


S12 ac^m4pivw 2 


Sj ac(,m4m^p2V 2 


s\2 acK,m.2m^p\vw 2 


sj ac( J m- i p4V 


■S6 «P3P5Z 


51 


s l2 ac 2 c 4 k 


S\2 ac^n^pivw 


2 2 

S12 acf > m4in^p2V w 


2 

s\2 ac 5112113 p\v w 


2 

J 12 aC(,m2P4VW 


s 6 ap 3 p 5 v 2 z 


52 


55 ac2Ct\k 


*i2 ac 3 n 4 p 1 


S\2 aC6«4«5P2 vvv 


•S12 ac 5 n 2 n i p 1 


s\2 ac^n^p4vw 


S12 a/?3P5VW 2 z 


53 


S13 a 2 C2C4\v 2 


2 

Sj a C~x,ilA,p\ 


•S12 ac 6 n 4 n 5 p 2 


2 

s-j a C5« 2 n3pi 


S12 ac 6 n 3 p 4 


■S12 ap 2 p 5 wz 


54 


4 

5 13 a C2C4W 


S13 a 3 c 3 ti4p\w 2 


sj a 3 c 6 n 4 n 5 p 2 


il3 a 3 C5«2 w 3PlW 2 


sj a 2 c 6 n 3 p 4 


■S12 aP3P5Z 


55 


*13 « 4 C2C4 


4 

S13 a C3/34P1W 


J13 a 4 c 6 n 4 n 5 p2W 2 


in a 5 c 5 n 2 n 3 piw 


s 13 a i C(,n^p4W 2 


s 6 a 3 p 3 p 5 z 


56 


s c 2 c 4 


s t3 a A c 3 n 4 pi 


J13 a 5 c 6 w 4 M 5 p2M' 


*13 a 5 c 5 n 2 n- i p x 


4 

i'13 a C(,n^p4\v 


■S13 a 3 P3P5W 2 


57 


•$0 C2C4 


*0 c 3Pl 


J13 a 5 c 6 n 4 n 5 p 2 




J13 a'C( s n 3 p4 


S13 a 3 P 3p 5 w 


58 


•So C2C4 


*0 c 3Pl 


•so t'6P2 


■SO c 5Pl 


S C 6 p 4 


S13 o 3 P3P5 
















59 


*0 C 2 C 4 


•SO c 3Pl 


■SO C(,P2 


J C'5Pl 


■SO C6P4 


■SO P3P5 



